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We study detailed classical-quantum correspondence for 
a cluster system of three spins with single-axis anisotropic 
exchange coupling. With autoregressive spectral estimation, 
we find oscillating terms in the quantum density of states 
caused by classical periodic orbits: in the slowly varying part 
of the density of states we see signs of nontrivial topology 
changes happening to the energy surface as the energy is var- 
ied. Also, we can explain the hierarchy of quantum energy lev- 
els near the ferromagnetic and antiferromagnetic states with 
EKB quantization to explain large structures and tunneling 
to explain small structures. 



I. INTRODUCTION 

When S is large, spin systems can be modeled by clas- 
sical and semiclassical techniques. Here we reserve "semi- 
classical" to mean not only that the technique works in 
the limit of large S (as the term is sometimes used) but 
that it implements the quantum-classical correspondence 
(relating classical trajectories to quantum-mechanical be- 
havior) . 

Spin systems (in particular S = 1/2) are often thought 
as the antithesis of the classical limit. Notwithstanding 
that, classical-quantum correspondence has been studied 
at large values of S in systems, such as an autonomous 
single spinid, kicked single spina, and autonomous twoti 
and threea spin systems. 

When the classical motion has a chaotic regime, for 
example, the dependence of level statistiosj-pn the regu- 
larity of classical motion has been studiedcllj. In regimes 
where the motion is predominantly regular, the pattern 
of quantum levels of a spin cluster can be understood 
with a combination of EBK (Einstein-Brillouin-Keller, 
also called Bohr-Sommerfeld) quantization and tunnel 
splitting (Sec. ^ is a such a study for the current system.) 
The latter sort of calculation has potential applications 
to some problems of current numerical or experimental 
interest. Numerical diagonalizations for extended spin 
systems (in ordered phases) on lattices of modest size (10 
to 36 spins) may be analyzed by treating the net spin of 
each sublattice as a single large spin and thereby reduc- 
ing the system to an autonomous cluster of a few spins; 
the clustering of low-lying eigenvalues can probe symme- 
try breakings that are obscured in a system ofrSuch size 
if only ground-state correlations are examined.^ Nonlin- 
ear self-localized modes in spin latticesQ, which typically 
span several sites, have to date been modeled classically, 
but seem well suited to semiclassical techniques. Another 



topic of recent experiments is the molecular mag netsB 
such as Mni2Ac and Fcg, which are more precisely mod- 
eled as clusters of several interacting spins rather than a 
single large spin; semiclassical analysis may pmvide an 
alternative to exact diagonalization techniquesO for the- 
oretical studies of such models. 

In this paper, we will study three aspects of the 
classical-correspondence of an autonomous cluster of 
three spins coupled by easy-plane exchange anisotropy, 
with the Hamiltonian 



H = 



(1) 



This model was introduced in Ref. ^, a study of level 
repulsion in regions of {E, a) space, where the classical 
dynamics is predominantly chaoticu. Eq. (|l|) has only 
two nontrivial degrees of freedom, since it conserves total 
angular momentum around the z-axis. As did Ref. ^ 
we consider only the case of Sf = 0. While studying 
classical mechanics we set |S| = 1; to compare quantum 
energy levels at different 5, we we divide energies by by 
S{S + 1) to normalize them. The classical maximum 
energy, E = 3, occurs at the ferromagnetic (FM) state 
- all three spins are coaligned in the equatorial (easy) 
plane. The classical ground state energy is E = —1.5, in 
the antiferromagnetic (AFM) state, in which the spins lie 
120° apart in the easy plane; there are two such states, 
differing by a reflection of the spins in a plane containing 
the z axis. Both the FM and AFM states, as well as all 
other states of the system, are continuously degenerate 
with respect to rotations around the z-axis. The classical 
dynamics follows from the fact that cos^; and 4>i, are 
conjugate, where 9i and (jji are the polar angles of the 
unit vector S^; then Hamilton's equations of motion say 



dcosOi/dt - 

d(l),/dt = ~h-'dHld(j), 



h-'dH/d4>i; 

-1 . 



(2) 



In the rest of this paper, we will first introduce the 
classical dynamics by surveying the fundamental periodic 
orbits of the three-spin cluster, determined by numeri- 
cal integration of the equations of motion (Sec |l^. The 
heart of the paper is Sec. IH: starting from the quantum 
density of states (DOS) obtained from numerical diago- 
nalization, we apply nonlinear spectral analysis to detect 
the oscillations in the quantum DOS caused by classical 
periodic orbits; to our knowledge, this is the first time 
the DOS has been related to specific orbits in a multi- 
spin system. Also, in Sec. IV we smooth the DOS and 



compare it to a lowest-order Thomas-Fermi approxima- 
tion counted by Monte Carlo integration of the classical 



1 



2 



energy surface; a flat interval is visible in the quantum 
DOS between two critical energies where the topology of 
the classical energy surface changes. Finally, in Sec. 0, 
we use a combination of EBK quantization and tunneling 
analysis to explain the clustering patterns of the quan- 
tum levels in our system. 



II. CLASSICAL PERIODIC ORBITS 

Our subsequent semiclassical analysis will depend on 
identification of all the fundamental orbits and their qual- 
itative changes as parameters are varied. Examining 
Poincare sections and searching along symmetry lines of 
the system, we found four families of fundamental pe- 
riodic orbits for the three-spin cluster. Figure |l| is an 
illustration of their motion, and Figure || gives classical 
energy-time curves. Orbits of types (a)-(c) are always 
at least threefold degenerate, since one spin is different 
from the other two; orbits of types (a)-(c) are also time- 
reversal invariant. Orbit (a), the counterbalanced orbit, 
exists when E > —1 (including the FM limit) and, in the 
range < a < 1 which we've studied, is always stable. 
Orbit (b), the unbalanced orbit, is unstable and exists 
when E < Ep, where 
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(3) 



Orbits of type (c), or stationary spin exist at all energies. 
Type (c) orbits are are unstable in the range, —1>E> 
3. Below E = —1 the stationary spin orbit bifurcates 
into two branches without breaking the symmetry of the 
ferromagnetic ground state. At 



(4) 



one branch vanishes and the other branch bifurcates into 
two orbits that are distorted spin waves of the two AFM 
ground states. (Below, in Sections IV and we will 
discuss topology changes of the entire energy surface.) 
Although they are not related by symmetry, all orbits of 
type (c) at a particular energy have the same period. 

Orbits of type (d), or three-phase orbits are named 
in analogy to three-phase AC electricity, as spin vectors 
move along distorted circles, 120° out of phase. The type 
(d) orbits break time-reversal symmetry and are hence 
at least twofold degenerate. A symmetry-breaking pitch- 
fork bifurcation of the (d) family occurs (for a = 0.5 
around E = —0.75) at which a single stable orbit, ap- 
proaching from high energy, bifurcates into an unstable 
and two stahle^recessing three-phase orbits without pe- 
riod doubling.llJ (Strictly speaking, the precessing three- 
phase orbits are not periodic orbits of the three-spin sys- 
tem, since after one "period" the spin configuration is 
not the same as before, but rather, all three spins are 
rotated by the same angle around the z-axis). The un- 
stable three-phase orbit disappears quickly as we lower 



energy, but the precessing three-phase orbits persist until 
E = —1.5, and become intermittently stable and unstable 
in a heavily chaotic regime near Ea, but regain stabil- 
ity before E —1.5: thus in the AFlVLJimit, orbits (c) 
are stable while orbits (b) are unstablcllJ More informa- 
tion on the classical mechanics of this system appears in 
Refs. |ll| and || 



III. ORBIT SPECTRUM ANALYSIS 

Gutzwiller's traejs formula, the central result of peri- 
odic orbit theory,t2l 

p{E)=ReJ2ME)eM^Sp{E)/h]+ptfiE), (5) 



decomposes the quantum DOS p{E) into a sum of oscil- 
lating terms contributed by classical orbits indexed by 
p, where Sp{E) is the classical action, and Ap{E) is a 
slowly jjizarying function of the period, stability and geo- 
metricLj properties of the orbit p) , plus the zeroth-order 
Thomas-Fermi term. 
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(6) 



This integral over phase space z is simply proportional 
to the area of the energy surface. We do not know of 
any mathematical derivation of (|^) in the case of a spin 
system. 

At a fixed H, the orbit spectrum is, as function of r, 
the power spectrum of p{E) inside the energy window, 
H-AH/2 <E < H+AH/2. (Figure|, explained below, 
is an example of an orbit spectrum.) Since the classical 
period Tp{E) = dSp{E)/dE, Eq. (|) implies that 0(iJ, r) 
is large if there exists a periodic orbit with energy H and 
period r. ThjC orbit spectrum can be estimated by Fourier 
transformJl3 



0{H,t) = 



Variants of Eq. 



H+AH/2 



H-AH/2 



p{E)e 



''dE 



(7) 



have been used to extract infor- 



mation ajpo'it classical periodic orbits from quantum 
spectraOlla Unfortunately, the resolution of the Fourier 
transform is limited by the uncertainty principle, 5E5t = 
h/2. 

Nonlinear spectral estimation techniques, however, can 
surpass the resolution of the Fourier transformJ13 One 
such technique, harmonic inversion, has been successfully 
applied to scaling systemaH3 - i.e., systems like billiards 
or Kepler systems in which the (classical and quantum) 
dynamics at one energy are identical to those at any other 
energy, after a rescaling of time and coordinate scales. In 
a scaling system, windowing is unnecessary because there 
are no bifurcations and the scaled periods of orbits are 
constant. In this section, we will apply nonlinear spectral 
estimation to our system (^, which is nonscalingEi 
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A. Diagonalization 

To get the quantum level spectrum, we wrote software 
to diagonalize arbitrary spin Hamiltonians polynomial in 
(Sf,S^,Sf), where i is an index running over arbitrary 
A'' spins of arbitrary (and often large) spin S. The pro- 
gram, written in Java, takes advantage of discrete trans- 
lational and parity symmetries by constructing a basis 
set in which the Hamiltonian is block diagonal, letting 
us diagonalize the blocks independently with an opti- 
mized version of LAPACK. Picturing the spins in a ring, 
the Hamiltonian Eq. is invariant to cyclic permuta- 
tions of the spins, so the eigenstates are states of definite 
wavenumbera k = 0, ±^ (matrix blocks for k = ±^ 
are identical by symmetry) . In the largest system we di- 
agonalized (three-spin cluster with S — 65) , the largest 
blocks contained N — 4620 states. 



B. Autoregressive approach to construct spectrum 

The input to an orbit spectrum calculation is the list of 
discrete eigenenergies with total 5*2 = 0; no other infor- 
mation on the eigenstates (e.g. the wavenumber quantum 
number) is necessary. This level spectrum is smoothed 
by convolving with a Gaussian (width 10^'^ for Figure |^) 
and discretely sampling over energy (with sample spacing 
(5 = 4.5 X 10-'*.) 

We estimate the power spectrum by the autoregressive 
(AR) method. AR models a discretely sampled input 
signal, Ui (in our case the density of states) with a process 
that attempts to predict yi from its previous values, 

N 

Vi ^^aiy.i-j + x^. (8) 

Here TV is a free parameter which determines how many 
spectral peaks that model can fit; Refs. |l^ and dis- 
cuss guidelines for choosing N. Fast algorithms exist to 
implement least-squares, i.e. to choose N coefficients 
to minimize (within-constraints) X^^fi '^^ these we used 
the Burg algorithmllj. 

To estimate the power spectrum, we discard the orig- 
inal Xi and model Xi with uncorrelated white noise. 
Thinking of Eq. (^) as a filter acting on Xi, the power 
spectrum of yi is computed from the transfer function of 
Eq. (|) and is 

Unlike the discrete Fourier transform, P(z^) can be eval- 
uated at any value of v. In our application, of course, 6 
has units of energy, so v (more exactly v/h) actually has 
units of time and is to be identified with t in (fzh. 



C. Orbit spectrum results and discussion 

Figure |^ shows the orbit spectrum of our system with 
S' = 65 and a = 0.5. ; it is displayed as a 500 x 390 array of 
pixels, colored light where 0{H, r) is large. Each horizon- 
tal row is the power spectrum in an energy window cen- 
tered at H; we stack rows of varying H vertically. With 
a window width 250 energy samples long {SH = 0.1125,) 
we fit = 150 coefficients in Eq. (||). To improve visual 
resolution, we let windows overlap and spaced the centers 
of successive windows 25 samples apart. 

Comparing Figure I and Fi gure ^ we see that our orbit 
spectrum detects the fundamental periodic orbits as well 
as multiple transversals of the orbits. Interestingly, we 
produced Figure ^ before we had identified most of the 
fundamental orbits; Figure ^ correctly predicted three 
out of four families of orbits. 

We believe that, given the same data, the AR method 
normally produces a far sharper spectrum. This is not 
surprising, since the Fourier analysis allows the possibil- 
ity of orbit-spectrum density at all r values, whereas AR 
takes advantage of our a priori knowledge that there are 
only a few fundamental periodic orbits and hence only 
a few peaks. We have compared the Fourier and AR 
versions of the spectrum in a few cases, but have not 
systematically tested them against each other. 

Unfortunately, the artifacts and limitations of the AR 
method are less understood than those of the Fourier 
transform. At high energies, the classical periods are 
nearly degenerate, so we expect closely spaced spectral 
peaks in the orbit spectrum. In this situation, the Burg 
algorithm vacillates between fitting one or two peaks 
causing the braiding between the (a) and (c) orbits (la- 
beled in Figure ^ in Figure |3[ Also, in the range 
— 1 < E < —1.3, where classical chaos is widespread, bi- 
furcations increase the number of contributing orbits so 
that we cannot interpret the orbit spectrum for r > 10. 

IV. AVERAGED DENSITY OF STATES 

The lowest-order Thomas Fermi approximation, Eq. 
(^) predicts that the area of the classical energy surface 
is proportional to the DOS. We verify this in Figure 
a comparison of the heavily smoothed quantum DOS to 
the area of the energy surface computed by Monte Carlo 
integration. 

An energy interval is visible in which the quantum 
DOS appears to be constant; we then verified that the 
classical DOS (which is more precise) is constant to our 
numerical precision; a similar interval was observed for 
all values of a. We identified this interval as {Ep,—1), 
where the endpoints are associated with changes in the 
topology of the energy surface as the energy varies. 

At energies below Ec (see Eq. (^)), the energy sur- 
face consists of two disconnected pieces, one surrounding 
each AFM ground state. The two parts coalesce as the 
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energy surface becomes multiply connected at Ec- For 
E < Ep, (see Eq. (^)) the anisotropic interaction con- 
fines the spins to a limited band of latitude away from 
the poles. At Ep it becomes possible for spins to pass 
over the poles. At E = —1, the holes that appeared in 
the energy surface at Ec close up. A discontinuity in the 
slope of the area of the energy surface occurs at energy 
Ec (not visible in Figure P; in the range Ep < E < —1 
the area of the energy surface (and hence the slowly vary- 
ing part of the DOS) seems to be constant as a function 
of energy. 

In the special isotropic (cr — 0) case, the flat interval 
is (—1.5, —1) and it can be analytically derived that the 
DOS is constant there. This is simplest for the smoothed 
quantum DOS, since for n = 1,2,... there are clusters 
of n energy levels with level spacing proportional to n. 
(A derivation also exists for the classical case, but is less 
direct.) We have no analytic results for general a. 

This flat interval is specific to our three-spin cluster, 
but we expect that the compactness of spin phase space 
will, generally, cause changes in the energy surface topol- 
ogy of spin systems that do not occur in traditionally 
studied particle systems. 

V. LEVEL CLUSTERING 

The quantum levels with total = show rich pat- 
terns of clustering, some of which arc visible on Figure |^. 
The levels that form clusters correspond to three differ- 
ent regimes of the classical dynamics in which the motion 
becomes nearly regular: (1) the FM limit (not visible in 
Figure ||; (2) the AFM limit (bottom edge of Figure Jp) 
and (3) the isotropic Hmit tr = (left edge of Figure 
Indeed, the levels form a hierarchy in as the clusters break 
up into subclusters. In this section, we first approxi- 
mately map the phase space from four coordinates to 
two coordinates - with the topology of a sphere. (Two of 
the original six coordinates are trivial, or decoupled, due 
to symmetry, as noted in Sec. ||. Then, using Einstein- 
Brillouin-Kramers (EBK) quantization some considera- 
tion of quantum tunneling, many features of the level 
hierarchy will be understood. 

A. Generic behavior: the polyad phase sphere 



In all three limiting regimes, the classical dynamics 
becomes trivial. For small deviations from the limit, the 
equations of motion can be linearized and one finds that 
the trajectory decomposes into a linear combination of 
two harmonic oscillators with degenerate frequency w, 
i.e., in a 1:1 resonance; the oscillators are coupled only 
by higher-order (=nonlinear) terms. 

There is a general prescription far understanding the 
classical dynamics in this situationo. Near the limit, the 



low-excited levels have approximate quantum numbers 
ni^2 such that the excitation energy AEi in oscillator i is 
hu!{ni + l/2). (In the FM limit, regime (1), this difference 
is actually measured downwards from the energy maxi- 
mum.) Clearly, the levels with a given total quantum 
number P = (rii -|- ?t.2 + 1) must have nearly degener- 
ate energies, and thus form a cluster of levels, which are 
split only by the effects (to be considered shortly) of the 
anharmonic perturbation— A level cluster arising in this 
fashion is called a polya(t3. 

To reduce the classical dynamics, make a canonical 
transformation to the variables <& and P = {Px,Py,Pz), 
where $ is the mean of the oscillators' phases and "ii^ is 
their phase difference, and 

Px = ^("-1 - n2), 
{Py,P,) = 2V(ni + l/2)(n2 + l/2)(cos*„sin*,), (10) 

Here $ is the fast coordinate, with trivial dynamics 
d^/dt — uj in the harmonic limit. The slow coordi- 
nates P follow a trajectory confined to the "polyad phase 
sphere" |P| = P, since AE = huP is conserved by the 
harmonic-order dynamics. The reduced dynamics on this 
sphere is properly a map P^ Pi+i, defined by (say) 
the Poincare section at '^x — O(mod 2tt). But dP/dt 
contains only higher powers of the components of P, so 
near the harmonic (small P) limit, |Pi+i — Pij-yanishes 
and the reduced dynamics becomes a flow. Lj At the 
limit in which it is a flow, an effective Hamiltonian / can. 
be defined so that the dynamics becomes integrable.Ej 
Applying EBK quantization to the reduced dynamics on 
the polyad phase sphere gives the splitting of levels within 
a polyad cluster. (Near the harmonic limit, the energy 
scale of / is small compared to the splitting between 
polyads.) 

In all three of our regimes, we believe this flass; has 
the topology shown schematically in Figure ^.lj Be- 
sides reflection symmetry about the "equator" , it also 
has a threefold rotation symmetry around the Pz axis, 
whichpCorrcsponds to the cyclic permutation of the three 
spins. Lj (Figure || is natural for the three-spin system 
because it is the simplest generic topology of the phase 
sphere with that threefold symmetry.) The reduced dy- 
namics has two symmetry-related fixed points at the 
"poles" Pz — ±P, which always correspond to motions 
of the three-phase sort like (d) on Figure |l|. There are 
also three stable and three unstable fixed points around 
the "equator". 

The K AM tori of the full dynamics correspond to orbits 
of the reduced dynamics. These orbits follow contours of 
the effective Hamiltonian / of the reduced dynamics (as 
in Figure ^). In view of the symmetries mentioned, 

I = aP^ + f3{P!^ - 3PxPy) + const (11) 

to leading order, where a,P, and the constant may de- 
pend on (7, S, and P. 
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The KAM tori surrounding the three-phase orbits rep- 
resented by the "poles" are twofold degenerate, while 
the tori in the stable resonant islands represented on the 
"equator" are threefold degenerate. Hence, the EBK con- 
struction produces degenerate subclusters containing two 
or three levels depending on the energy range within the 
polyad cluster. 

The fraction of levels in one or the other kind of sub- 
cluster is proportional to the spherical areas on the corre- 
sponding side of the separatrix, which passes through the 
unstable points in Figure ^. These areas in turn depend 
on the ratio of the first to the second term in Eq. (11), 
i.e. OiP^ I fiP^ . Evidently, as one moves away from the 
harmonic limit to higher values of P, one universally ex- 
pects to have a larger and larger fraction of threefold 
subclusters. 

Given the numerical values of energy levels in a polyad, 
we can estimate the terms of Eq. (nTl) in the following 
fashion: (i) the energy difference between the highest 
and lowest 3-fold subcluster is the difference between the 
stable and unstable orbits on the equator, which is 2/3P^ 
according to (pT|); (ii) the mean of the highest and low- 
est 3-fold subcluster would be the energy all around the 
equator if /3 were to vanish; the difference between this 
energy and that of the farthest 2-fold subcluster in the 
polyad is aP^ according to ([Til). 

Furthermore, tunneling between nearby tori creates 
fine structure splitting inside the sub-clusters. The slow 
part of the dynamics on the polyad phase sphere, is iden- 
tical to that of a single semiclassical spin with ( |lTl) as its 
effective Hamiltonian, so the effective Lagrangian is es- 
sentially the same, too. Then different tunneling paths 
connecting the same two quantized orbits must differ in 
phase by a topological term, with a familiar form pro- 
portional to the-|(real part of the) spherical area between 
the two paths. Ea 



k — ±1 and every "1" lev el has k 
ber k was defined in Sec. 
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unstable points on the "equator" are identified respec- 
tively with the orbits (a) and (c) of Figure |^. Classically, 
in this regime, the three-phase orbit is the fundamen- 
tal orbit with lowest frequency io^-phase] thus the cor- 
responding levels in successive polyads have a somewhat 
smaller spacing huj^-phase than other levels, and they 
end up at the top of each polyad. (Remember excitation 
energy is measured downwards from the FM limit.) In- 
deed, we observe that the high-energy end of each polyad 
consists of twofold subclusters and the low-energy end 
consists of threefold subclusters. 

We see a pattern of fine structure (presumably tunnel 
splittings^ r-which is just like the pattern in the four-spin 
problemB'CJ Namely, throughout each polyad the degen- 
eracies of successive levels follow the pattern (2,1,1,2) and 
repeat. (Here - as also for regime 3 - every "2" level has 

= 0, where wavenum- 
Numerical data show 
that (independent of S) the pattern (starting from the 
lowest energy) begins (2112...) for even P, but for odd 
P it begins (1221...). 

In the energy range of twofold subclusters, the levels 
are grouped as (2) (11) (2), i.e. one tunnel-split subclus- 
ter between two unsplit subclusters(and repeat); in the 
threefold subcluster regime, the grouping is (21)(12), so 
that each subcluster gets tunnel-split into a pair and a 
single level, but the sense of the splitting alternates from 
one subcluster to the next. 

An analysis of u = 0.4, 5 = 30 showed that the fraction 
of threefold subclusters indeed grows from around 0.3 for 
small P to nearly 0.5 at P w 40. Furthermore, when aP^ 
and l3P^ were estimated by the method described near 
the end of Subscc. V A, they indeed scaled as P^ and P^ 
respectively. 



2. Antiferromagnetic limit 



B. Results 

Here we summarize some observations made by exam- 
ination of polyads in the three regimes, for a few combi- 
nations of S and a. 



1. Ferromagnetic limit 

This regime is the best-behaved in that regular behav- 
ior persists for a wide range of energies. The ferromag- 
netic state, an energy maximum, is a fixed point of the 
dynamics; around it are "spin-wave" excitations (viewing 
our system as the 3-site case of a one dimensional ferro- 
magnet). These are the two oscillators from which the 
polyad is constructed. Thus, the "pole" points in Fig- 
ure H correspond to "spin waves" propagating clockwise 
or counterclockwise around the ring of three spins, an ex- 
ample of the "three-phase" type of orbit. The stable and 



This regime occurs at _B < Ec{(t), where Ec{a) is given 
by (^ . That means the classical energy surface is divided 
into two disconnected pieces, related by a mirror reflec- 
tion of all three spins in any plane normal to the easy 
plane. Analogous to regime one, two degenerate antifer- 
romagnetic "spin waves" exist around either energy min- 
imum, and the polyad states are built from the levels of 
these two os cilla tors. Thus the clustering hierarchy out- 
lined in Sec. V A - polyads clusters, EBK-quantization of 



/, and tunneling over barriers of / on the polyad phase 
sphere - is repeated within each disconnected piece, lead- 
ing to a prediction that all levels should be twofold de- 
generate. 

Consequently, on the level diagram (Figure ^), there 
should be half the apparent level density below the line 
E = Ec{iy) as above it. Indeed, a striking qualitative 
change in the apparent level crossing behavior is visible 
at that line (shown dashed in the figure). 

Actually, tunneling is possible between the discon- 
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nected pieces of the energy surface and may split these 
degenerate pairs. In fact this hyperfine sphtting hap- 
pens to 1/3 of the pairs, again following the (2112) pat- 
tern within a given polyad. This (2112) pattern starts to 
break up as the energy moves away from the AFM limit; 
even for large S (30 or 65), this breakup happens already 
around the polyad with P = 10, so it is much harder 
than in the FM case to ascertain the asymptotic pattern 
of subclustering. We conjecture that the breakup may 
happen near the energies where, classically, the stable 
periodic orbits bifurcate and a small bit of phase space 
goes chaotic. 

The barrier for tunneling between the disconnected en- 
ergy surfaces has the energy scale of the bare Hamilto- 
nian, which is much larger (at least, for small P) than 
the scale of effective Hamiltonian / which provides the 
barrier for tunneling among the states in a subcluster. 
Hence, the hyperfine splittings are tiny compar ed to the 
fine sphttings discussed at the end of Subsection V A . To 
analyze numerical results, we replace a degenerate level 
pair by one level and a hyperfine-split pair by the mean 
level, and treat the result as the levels from one of the two 
disconnected polyad phase spheres, neglecting tunneling 
to the other one. 

Then in the AFM limit, the "pole" points in Figure ^ 
again correspond to spin waves propagating around the 
ring, while the stable and unstable points on the equa- 
tor are (c) and (b) on Figure |l]. The three-phaSje-|Or- 
bit is the highest frequency orbit in the AFM limit £3 so 
again the twofold and threefold subclusters should occur 
at the high and low energy ends of each polyad cluster. 
What we observe, however, is that all the subclusters are 
twofold, except the lowest one is often threefold. 



3. Isotropic limit 

This regime will includes only Stot < S i.e. E < ~1 
- the same regime in which the flat DOS was observed 
(Sec. IV). Above the critical value E = —1, the levels 
behave as in the "FM limit" described above. 

At cr = 0, it is well-known that the quantum Hamilto- 
nian reduces to ^[S^^^ — 3S{S + 1)]. Thus each level has 
degeneracy P = 2S'tot + 1- (That is the number of ways 
three spins S may be added to make total spin Stot , and 
each such multiplet has one state with St^t = 0.) Whea, 
a is small, these levels split and will be called a polyad.E2l 

Classically, at ct = the spins simply precess rigidly 
around the total spin vector. These are harmonic mo- 
tions of four coordinates; hence the polyad phase sphere 
can be constructed by (p^). From the threefold symme- 
try, there should again be three orbit types as represented 
generically by Figure |^ and Eq. (^Tj). For example, an 
umbrella-like configuration in which the three spin direc- 
tions are equally tilted out of their plane corresponds to a 
three-phase type orbit, with two cases depending on the 
handedness of the arrangement. A configuration where 



one spin is parallel/ antiparallel to the net moment, and 
the other two spins offset symmetrically from it), follows 
one of the threefold degenerate orbits. 

Numerically, the level behavior in the near-isotropic 
limit is similar to the near-FM limit. The fine structure 
degeneracies are a repeat of the (2112) pattern as in the 
other regimes; the lowest levels of any polyad always be- 
gin with (1221). The fraction of threefold subclusters is 
large here and, as expected, grows with S, (from 0.5 to 
0.7 in the case S = 15). However, the energy scales of 
aP^ and PP^ behave numerically as aP'^ and aP^ . What 
is different about the isotropic limit is that the precession 
frequency - hence the oscillator frequency a; ~ is not a 
constant, but is proportional to Stot- Since perturbation 
techniquesLj give formulas for / with inverse powers of 
w, it is plausible that a and /3 in (jl^) include factors of 
P^^ here, which were absent in the other two regimes. 



VI. CONCLUSION AND SUMMARY 

To summarize, by using detailed kn|©wledge of the clas- 
sical mechanics of a three spin clusteiO, we have studied 
the semiclassical limit of spin in three ways. First, using 
autoregressive spectral analysis, we identified the oscil- 
lating contributions that the fundamental orbits of the 
cluster make to the density of states, in fact, we detected 
the quantum signature of the orbits before discovering 
them. Secondly, we verified that the quantum DOS is 
proportional to the area of the energy surface; we also 
observed kinks in the smoothed quantum DOS, which 
are the quantum manifestation of topology changes of 
the classical energy surface; such topology changes, we 
expect, are more common in spin systems than parti- 
cle phase space, since even a single spin has a nontriv- 
ial topology. Finally, we have identified three regimes of 
near-regular behavior in which the levels are clustered ac- 
cording to a four-level hierarchy, and we explained many 
features qualitatively in terms of a reduced, one degree- 
of-freedom system. This system appears promising for 
two extensions analgous to Ref. |5|: tunnel amplitudes 
(and their topological phases) could be computed more 
explicitly; also, the low-energy levels from exact diag- 
onalization of a finite piece of the anisotropic-exchange 
antiferromagnet on the triangular lattice could probably 
be mapped to three large spins and analyzed in the fash- 
ion sketched above in Sec. v[ 
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FIG. 3. Orbit spectrum for S — 65, a = 0.5. The horizontal 
FIG. 1. Fundamental periodic orbits of the tSree-spin clus- axis is classical period, and the vertical axis is energy. Peaks 
ter. (a) Counterbalanced, (b) unbalanced orbits, (c) station- ^he orbit spectrum are white and valleys are black, 

ary spin, and (d) three-phase. 

80.0 r 




period _i.o -0.5 0.0 



FIG. 2. Energy-period curve of the three spin system energy 

with a = 0.5. The curves (a), (b), (c) and (d) are four fami- FIG. 4. The density of states of a 5 = 40 three spin cluster 

lies of periodic orbits . Below the obvious bifurcation around with a = 0.1, smoothed with a Gaussian exp[—{E/S)^] with 

E = —0.75, orbit (d) is not literally a periodic orbit of the 5 = 0.01. The Thomas-Fermi density of states is very flat in 

three spins, but only in a reduced two degree of freedom sys- the range Ep < E < —1, with Ep = —1.43. 
tem (wherein one identifies states related by a rotation of all 
three spins about the z axis. 



-1.50 



FIG. 5. Energy levels as a function of a for an S = 10 
three spin cluster. The dashed line is the coalescence energy 
Ec{<j) defined by equation Eq. (0). 




FIG. 6. Topology of orbits for a generic phase sphere 
with 3-fold symmetry. The phase space is decomposed into 
one "fast" action variable, which is transverse to the 3-space 
of the sphere, and three "slow" variables; of these, the en- 
ergy is a function of radius and is conserved. Contours are 
shown of an effective Hamiltonian /, with the direction of 
"slow" motion indicated by arrows. Tunneling takes place 
between symmetry-related orbits across the saddle points of 
/, as indicated by the dashed lines. The stationary points of 
this two-dimensional flow labeled (b), (c), and (d) correspond 
near the antiferromagnetic ground state to the periodic or- 
bits of the full system labeled the same way in Figure ^; the 
labeling would differ somewhat for t he near-ferromagnetic or 
near- isotropic regimes (see Sec. VB). 



